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ABSTRACT 

Context. The Hundred-Thousand-Proper-Motion (HTPM) project will determine the proper motions of ~1 13, 500 stars using a ~23- 
year baseline. The proper motions will be based on space-based measurements exclusively, with the Hipparcos data, with epoch 
1991.25, as first epoch and with the first intermediate-release Gaia astrometry, with epoch ~2014.5, as second epoch. The expected 
HTPM proper- motion standard errors are 30-190 yuas yl -1 , depending on stellar magnitude. 

Aims. Depending on the astrometric characteristics of an object, in particular its distance and velocity, its radial velocity can have a 
significant impact on the determination of its proper motion. The impact of this perspective acceleration is largest for fast-moving, 
nearby stars. Our goal is to determine, for each star in the Hipparcos catalogue, the radial- velocity standard error that is required to 
guarantee a negligible contribution of perspective acceleration to the HTPM proper-motion precision. 

Methods. We employ two evaluation criteria, both based on Monte-Carlo simulations, with which we determine which stars need to 
be spectroscopically (re-)measured. Both criteria take the Hipparcos measurement errors into account. The first criterion, the Gaussian 
criterion, is applicable to nearby stars. For distant stars, this criterion works but returns overly pessimistic results. We therefore use a 
second criterion, the robust criterion, which is equivalent to the Gaussian criterion for nearby stars but avoids biases for distant stars 
and/or objects without literature radial velocity. The robust criterion is hence our prefered choice for all stars, regardless of distance. 
Results. For each star in the Hipparcos catalogue, we determine the confidence level with which the available radial velocity and its 
standard error, taken from the XHIP compilation catalogue, are acceptable. We find that for 97 stars, the radial velocities available in 
the literature are insufficiently precise for a 68.27% confidence level. If requiring this level to be 95.45%, or even 99.73%, the number 
of stars increases to 247 or 382, respectively. We also identify 109 stars for which radial velocities are currently unknown yet need to 
be acquired to meet the 68.27% confidence level. For higher confidence levels (95.45% or 99.73%), the number of such stars increases 
to 1, 071 or 6, 180, respectively. 

Conclusions. To satisfy the radial-velocity requirements coming from our study will be a daunting task consuming a significant 
amount of spectroscopic telescope time. The required radial-velocity measurement precisions vary from source to source. Typically, 
they are modest, below 25 km s -1 , but they can be as stringent as 0.04 km s _I for individual objects like Barnard's star. Fortunately, 
the follow-up spectroscopy is not time-critical since the HTPM proper motions can be corrected a posteriori once (improved) radial 
velocities become available. 

Key words. Techniques: radial velocities, Astronomical databases: miscellaneous, Catalogs, Astrometry, Parallaxes, Proper motions 



1. Introduction 



Gaia (e.g., |Perryman et"aL"l|2001 |Lindegren et al. 2008) is the 
upcoming astrometry mission of the European Space Agency 
(ESA), following up on the success of the Hipparcos mission 
pSA] |1997a[ |Perryman efaL] [T997} |Perryman| [2009] i. Gaia's 
science objective is to unravel the kinematical, dynamical, and 
chemical structure and evolution of our galaxy, the Milky Way 
(e.g., |G6mez et al.| 



addition, Gaia's data will revo- 
lutionise many other areas of (astro )physics, e.g., stellar struc- 
ture and evolution, stellar variability, double and multiple stars, 
solar-system bodies, fundamental physics, and exo-planets (e.g., 

~ [2008) |Mignard & Klioner] [jolTj 



Pourbabq 2008 ; Tanga et al. 



Eyer et al.||201 i[|Sozzetti[|20i l||Mouret||201 1| >. During its five 



year lifetime, Gaia will survey the full sky and repeatedly ob- 
serve the brightest 1,000 million objects, down to 20 th magni- 
tude (e.g., |de Bruijne et al.||2010| . Gaia's science data comprises 
absolute astrometry, broad-band photometry, and low-resolution 
spectro-photometry. Medium-resolution spectroscopic data will 
be obtained for the brightest 150 million sources, down to 17 th 



magnitude. The final Gaia catalogue, due in ~2021, will con- 
tain astrometry (positions, parallaxes, and proper motions) with 
standard errors less than 10 micro-arcsecond (yuas, yuas yr~' for 
proper motions) for stars brighter than 12 th magnitude , 25 ;uas for 
stars at 15 th magnitude, and 300 yuas at magnitude 20 ( de Bruijne| 
|2012J ). Milli-magnitude-precision photometry ( |Jordi et aT 2010) 
allows to get a handle on effective temperature, surface grav- 
ity, metallicity, and reddening of all stars (Bail er-Jones| |2010| l. 
The spectroscopic data allows the determination of radial veloc- 
ities with errors of 1 km s" 1 at the bright end and 15 km s 
at magnitude 17 ( jWilkinson et al.j |2005[ |Katz et al.[ |2011| i as 
well as astrophysical diagnostics such as effective temperature 
and metallicity for the brightest few million objects ( |Kordo patis 
|et al.| |2011| l. Clearly, these performances will only be reached 
with a total of five years of collected data and only after careful 
calibration. 

Intermediate releases of the data - obviously with lower 
quality and/or reduced contents compared to the final cat- 
alogue - are planned, the first one around two years after 
launch, which is currently foreseen for the second half of 
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2013. The Hun dred-Thousand-Proper-Motion (HTPM) project 
(Mignard 2009), conceived and led by Francois Mignard at the 
Observatoire de la Cote d'Azur, is part of the first intermedi- 
ate release. Its goal is to determine the absolute proper motions 
of the ~1 13, 500 brightest stars in the sky using Hipparcos as- 
trometry for the first epoch and early Gaia astrometry for the 
second. Clearly, the HTPM catalogue will have a limited life- 
time since it will be superseded by the final Gaia catalogue 
in ~2021. Nevertheless, the HTPM is a scientifically interest- 
ing as well as unique catalogue: the ~23-year temporal base- 
line, with a mean Hipparcos epoch of 1991.25 and a mean Gaia 
epoch around 2014.5, allows a significant improvement of the 
Hipparcos proper motions, which have typical precisions at the 
level of 1 milli-arcsecond yr~' (mas yr -1 ): the expected HTPM 
proper-motion standard errorsH are 40-190 pas yr~' for the 
proper motion in right ascension fi a * and 30-150 pas yr~' for the 
proper motion in declination fig, primarily depending on magni- 
tude (we use the common Hipparcos notation a* = a ■ cos 5; 



ESA 1997a Section 1.2.5). A clear advantage of combining as- 



trometric data from the Hipparcos and Gaia missions is that the 
associated proper motions will be, by construction and IAU res- 
olution, in the system of the International Celestial Reference 
System (ICRS), i.e., the proper motions will be absolute rather 
than relative. In this light, it is important to realise that mas- 
sive, modern-day proper-motion catalogues, such as UCAC-3 
(Zacharias et al. |2010| l, often contain relative proper motions 
only and that they can suffer from substantial, regional, system- 
atic distortions in their prop er-motion systems, up to levels of 



2008 



2010 



Liu et al. 



10 mas yr 1 or more (e.g., |Roser et al. 
[20TT) . 

It is a well-known geometrical feature, for instance already 
described by Seeliger in 1900 that for fast-moving, nearby stars, 



it is essential to know the radial velocity for a precise mea- 
surement and determination of proper motion. In fact, this so- 
called secular or perspective acceleration on the sky was taken 
into account in the determination of the Hipparcos proper mo- 
tions for 21 stars ( |ESA| |1997a| Section 1.2.8) and the same 
will be done for Gaia, albeit for a larger sample of nearby 
stars. Clearly, the inverse relationship also holds: with a pre- 
cise proper motion available, a so-called astrometric radial ve- 
locity can be determined, independent of the spectroscopically 



measured quantity (see Lindegren & Dravins 2003 for a pre- 
cise definition and meaning of [astrometric] radial velocity). 
With this method, Dravins et al.| ( |1999| determinecj^] the astro- 
metric radial velocities for 17 stars, from Hipparcos proper mo- 
tions combined with Astrographic Catalogue positions at earlier 

1 The HTPM proper motions will be limited in precision by the 
Hipparcos parallax uncertainties, which are typically ~ 1 mas (the typ- 
ical HTPM proper-motion standard error is hence 1 mas/23 yr « 
40 yuas yr -1 ). The first intermediate-release Gaia catalogue is based on 
just ~12 months of data, which is generally insufficient to unambigu- 
ously lift the degeneracy between proper motion and parallax for all 
stars. The underlying astrometric global iterative solution (Lindegren 
et al. 2012) will hence be based on a two- rather than five-parameter 
source model, fitting for position (a, S) at the mean Gaia epoch only. 
The Hipparcos parallax is hence needed to correct the Gaia transit ob- 
servations for the parallactic effect allowing to transform apparent di- 
rections into barycentric positions. 

2 These authors also describe two other methods to derive astrometric 
radial velocities, namely by measuring changing annual parallax or by 
measuring changing angular extent of a moving group of stars (Madscn 
|et al.| 



paral 



et al. 



2002). The latter method also provides, as a bonus, improved 
axes to moving-group members (e.g., de Bruijne , 1999 , de Bruijne 



epochs. Although |Dravins et al.| ([T999 ) reached relatively mod- 
est astrometric-radial-velocity precisions, typically a few tens of 
km s _1 , their results are interesting since they provide direct and 
independent constraints on various physical phenomena affect- 
ing spectroscopic radial velocities, for instance gravitational red- 
shifts, stellar rotation, convection, and pulsation. In our study, 
however, we approach (astrometric) radial velocities from the 
other direction since our interest is to determine accurate HTPM 
proper motions which are not biased by unmodeled perspective 
effects. In other words: we aim to establish for which stars in 
the forthcoming HTPM catalogue the currently available (spec- 
troscopic) radial velocity and associated standard error are suf- 
ficient to guarantee, with a certain confidence level, a negligible 
perspective-acceleration-induced error in the HTPM proper mo- 
tion. For stars without a literature value of the radial velocity, we 
establish whether - and, if yes, with what standard error - a ra- 
dial velocity needs to be acquired prior to the construction of the 
HTPM catalogue. Section [2] describes the available astrometric 
and spectroscopic data. The propagation model of star positions 
is outlined in Section[3] We investigate the influence of the radial 
velocity on HTPM proper motions in Section[4]and develop two 
evaluation criteria in Section [5] We employ these in Section [6] 
We discuss our results in Section|7]and give our final conclusions 
in Section|8] 



2. The XHIP catalogue 

As source for the Hipparcos astrometry and literature radial 
velocities, we used the extended Hipparcos compilation cata- 



logue (CDS catalogue V/137), also known as XHIP (Anderson 



& Francis 2012| ). This catalogue complements the 117,955 en- 
tries with astrometry in the Hipparcos catalogue with a set of 
116,096 spectral classifications, 46,392 radial velocities, and 
18, 549 iron abundances from various literature sources. 



2.1. Astrometry 

The starting point for the XHIP compilation was the new reduc- 
tion of the Hipparcos data ( |van Leeuwen||200"7] [2008 , CDS cata- 
logue 1/3 11), also known as HIP-2. Realising that stars with mul- 
tiple components were solved individually, rather than as sys- 



tems, by |van Leeuwen for the sake of expediency, Ander son &| 
|Francis| reverted to the original HIP-1 astrometry (ESA| |1997a|b| 
CDS catalogue 1/239) in those cases where multiplicity is in- 
dicated and where the formal parallax standard error in HIP-2 
exceeds that in HIP-1. This applies to 1,922 entries. In addi- 
tion, Anderson & Francisl included the Tycho-2 catalogue (H0g 



2001) . 



|et aL| |2000a|b| CDS catalogue 1/259) in their XHIP proper- 
motion data. In the absence of Tycho-2 proper motions, HIP- 
2 proper motions were forcibly used. When multiplicity is in- 
dicated, Hipparcos proper motions were replaced by Tycho-2 
values in those cases where the latter are more precise. When 
multiplicity is not indicated, Tycho-2 proper motions replaced 
Hipparcos values if the associated standard errors exceed the 
Tycho-2 standard errors by a factor three or more. In all other 
cases, a mean HIP-2 - Tycho-2 proper motion was constructed 
and used, weighted by the inverse squared standard errors; this 
applies to 92, 269 entries. 

2.2. Radial velocities 

The XHIP catalogue contains radial velocities for 46, 392 of 
the 1 17, 955 entries, carefully compiled by |Anderson & Francisj 
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from 47 literature sources. The vast majority of measurements 
have formal measurement precisions, i.e., radial-velocity stan- 
dard errors (1,753 measurements lack standard errors; see 
Section [6!2| . In addition, all radial velocities have a quality flag: 

- An 'A' rating (35,932 entries) indicates that the standard er- 
rors are generally reliable. 

- A 'B' rating (4,239 entries) indicates potential, small, un- 
corrected, systematic errors. 

- A 'C rating (3,465 entries) indicates larger systematic er- 
rors, while not excluding suitability for population analyses. 

- A 'D' rating (2,756 entries) indicates serious problems, 
meaning that these stars may not be suitable for statistical 
analyses. A 'D' rating is assigned whenever: 

1. the radial-velocity standard error is not available, 

2. the star is an un(re)solved binary, 

3. the star is a Wolf-Rayet star or a white dwarf that is not 
a component of a resolved binary, or 

4. different measurements yield inconsistent results. 

The majority of stars in the XHIP catalogue (71, 563 entries 
to be precise) have no measured radial velocity. All we can rea- 
sonably assume for these stars is that their radial-velocity dis- 
tribution is statistically identical to the radial-velocity distribu- 
tion of the entries with known radial velocities. Figure [TJ shows 
this distribution. It is fairly well represented by a normal distri- 
bution with a mean fi = -2.21 km s and standard deviation 
cr = 22.44 km s _1 (the median is -2.00 km s _1 ). The observed 
distribution has low-amplitude, broad wings as well as a small 
number of real 'outliers', with heliocentric radial velocities up to 
plus-or-minus several hundred km s" 1 . A small fraction of these 
stars are early-type runawaj|^]stars (Hoogerwerf et al. 200 l| l but 
the majority represent nearby stars in the (non-rotating) halo of 
our galaxy. The bulk of the stars, those in the main peak, are 
(thin-)disc stars, co-rotating with the Sun around the galactic 
centre. In theory, the main peak can be understood, and also be 
modeled in detail and hence be used to statistically predict the 
radial velocities for objects without literature values, as a com- 
bination of the reflex of the solar motion with respect to the lo- 
cal standard of rest ( Delhaye 1965 Schonrich et al. 2010 1, the 

I1927| l, 



effect of differential galactic rotation (Oort 



dom motion of stars (ISchwarzschild 1907 



and the ran- 
In practice, how- 
ever, such a modeling effort would be massive, touching on a 
wide variety of (sometimes poorly understood) issues such as 
the asymmetric drift, the tilt and vertex deviation of the veloc- 
ity ellipsoid, mixing and heating of stars as function of age, the 
height of the sun above the galactic plane (Joshi, 2007), the dy- 
namical coupling of the local kinematics to the galactic bar and 
spiral arms ( Antoja et al. |201 l| l, large-scale deviations of the lo- 
cal velocity field caused by the Gould Belt (|Elias et al. ] [2006), 
migration of stars in the disc ( Schonrich & Binney , 2009), etc. To 
model these effects, and hence be able to predict a more refined 
radial velocity for any star as function of its galactic coordinates, 
distance, and age, is clearly beyond the scope of this paper. We 
will come back to this issue in SectionROl 



3. Propagation model 

3.1. The full model 

Let us denote the celestial position of a star at time to (in years) 
in equatorial coordinates in radians by (ao, 6q), its distance in 
parsec by r = 1000 • et , with the parallax zzr in mas, and its 




-200 -100 100 

Radial velocity [km s" 1 ] 



200 



Fig. 1. Distribution of all 46, 392 radial velocities contained in 
the XHIP catalogue. The smooth, red curve fits the histogram 
with a Gaussian normal distribution. The best-fit mean and stan- 
dard deviation are /i = -2.21 km s _1 and cr = 22 .44 km s , 
respectively. 

proper-motion components in equatorial coordinates in mas yr~' 
by (ji a ,fis)- The three-dimensional position of a star in Cartesian 
equatorial coordinates at time to is then given by: 

xq — r ■ cos a>o ■ cos So; 
yo = r ■ sin ao ■ cos do; 

zo = r ■ sin5 - (1) 

With v r the star's radial velocity in km s _1 , it is customary to 
define the 'radial proper motion' as = v r • r _1 . The linear 
velocity in pc yr~' is then given by: 

v a = y"ff • B ■ Ay ■ cr -1 ; 
Vg = [is • B ■ Ay ■ tu^ 1 ; 

V r = fl r ■ B, (2) 

where Ay = 4.740,470,446 km yr s 
and the factor B 
changes km s 



: A p -A- Z l 



is the astronomical unit 
1 .022, 712, 169 10 6 pc s knT 1 yr 1 



B -A v = 4.848,136,811 



1 to pc yr 1 ( |ESA| 1 1997a Table 1.2.2) so that 
lO^pc AU" 



Transforming these 



equations into Cartesian coordinates leads to: 

v x = — v a sinaro — v^cosao sin^o + v r cosao cos(5o; 
v a cos ao - vg sin ao sin do + v r sin ao cos <5o; 

\>s cosdo + V r sin<5rj. 



Vy 



(3) 



Since the motion of stars, or the barycentre of multiple systems, 
is to near-perfect approximation rectilinear over time scales of 
a few decades, the position of a star at time t\, after a time t - 
1 1 - to, now simply follows by applying the propagation model: 

x(t) = x + v x -t; 
y(t) = yo + v y ■ f; 

z(t) = zo + v z ■ t. (4) 

Transforming this back into equatorial coordinates returns the 
celestial position (a{t), 6{t)) of the star at time f: 



The Hipparcos catalogue does not contain hyper- velocity stars. 



a(t) = arctan 



5(t) = arctan 



x(t) 



z(t) 



V*(o 2 +y(0 2 



(5) 
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So, in summary, it is straightforward to compute the future po- 
sition of a star on the sky once the initial celestial coordinates, 
the proper motion, parallax, radial velocity, and time interval are 
given. This, however, is what nature and Gaia will do for us: the 
initial celestial coordinates correspond to the Hipparcos epoch 
(1991.25) and the final coordinates (a(t),6(t)), with t w 23 yr, 
will come from the first Gaia astrometry. The HTPM project will 
determine the proper- motion components (ii a ,[i$) from known 
initial and final celestial coordinates for given time interval, 
parallax, and radial velocity. Unfortunately, it is not possible 
to express the proper-motion components in a closed (analyti- 
cal) form as function of (ao,do), (a(t),S(t)), m, v r , and t since 
the underlying set of equations is coupled. The derivation of 
the proper-motion components hence requires a numerical so- 
lution. We implemented this solution using Newton-Raphson 
iteration and refer to this solution as the 'full model'. This 
model, however, is relatively slow for practical implementation. 
We hence decided to also implement a 'truncated model' with 



analytical terms up to and including t (Section 3.2 1, which is 



about ten times faster and sufficiently precise for our application 
(SectionfO). 



3.2. The truncated model 

Mignard (2009) shows that the full model (Section [3. l| l can be 
truncated up to and including third-order terms in time t without 
significant loss in accuracy (Section 3.3 I. Equations |6]) and (|7]i 
below show the forward propagation for right ascension and dec- 
lination, respectively. By forward propagation, we mean com- 
puting the positional displacements Aa and A6 of a star for 
given proper motion, parallax, radial velocity, and time interval 
t = h — to'- 

Aar* = Aar cos 6q = (a(t) - a?o) cos 6q = 

\Ma] ■ t 1 + 
- [yu,.yu ff - tanSofia/ug] ■ t 2 + 



l4l* a - 2 tan SofJ-rVaVs + tan 2 6 (J. a (4 



3 cos 2 6[) 



f 3 + 



0(f); (6) 



AS = 6(t) -6 = 
M ■ t 1 + 



tan 6o 2 



/4ns + tan 5ofi r nl - 



2 cos 2 do 



t 2 + 



r 3 + 



0(t 4 ). 



(7) 



Equations ([8]l and (|9]l below show the backward solution for the 
proper motion in right ascension and declination, respectively. 
By backward solution, we mean computing the proper-motion 
components (fi a ,fi$) from the initial and final celestial positions 
(ao,<5o) and (a(t),6(t)), for given parallax, radial velocity, and 
time interval t — t\ - to: 

fJ-a-t = 

Aa* + Aa*f£ r t + 
- tan SoAa* A6 - tan 6oAa*A6fi r t + 
3 cos 2 do - 1 



t*s ■ t = 

Ad + A6/u r t + 

1 ,1 , 

+ - tm6o(Aa*Y + - tanSo(Aa*r/j. r t + 

2cos 2 <5 -l,. , v2 .. A<5 3 
2 cos z do 3 



(9) 



It is straightforward to insert Equations (|6]l and |7]i into 
Equations ([8]l and (|9|l to demonstrate that only terms of order 
f 4 and higher are left. 



3.3. Accuracy of the truncated model 

To quantify that the truncation of the full model up to and includ- 
ing third-order terms in time is sufficient for the HTPM applica- 
tion, Table [T] shows the errors in derived proper motions over 
an interval of 25 years induced by the truncation of the model 
when using the approximated Equations |6]l-((9]l for an 'extreme' 
star (i.e., nearby, fast-moving and hence sensitive to perspective- 
acceleration effects) as function of declination. Four cases have 
been considered. First, Equations |6]l-((7]) up to and including 
first-order terms in time were used for the forward propagation 
and Equations ([8])-((9]) were used for the backward solution. This 
is indicated by the heading 0(t l ). The difference between the 
proper motion used as input and the proper motion derived from 
Equations ([8]l-((9]) is listed in the table and can reach several 
mas yr~' close to the celestial poles. The second case ('(9(f 2 )') is 
similar to the first case but includes also second-order terms in 
time for the forward propagation. The proper-motion errors are 
now much reduced, by about three orders of magnitude, but can 
still reach 10 yuas yr~', which is significant given the predicted 
HTPM standard errors (30-190 //as yr -1 , depending on magni- 
tude). The third case ('(9(f 3 )') also includes third-order terms 
in time. The proper-motion errors are now negligible, reaching 
only up to 10 nano-arcsecond yr~'. The fourth case uses the full 
model for the forward propagation and the Newton-Raphson it- 
eration for the backward solution and recovers the input proper 
motions with sub-nano-arcsecond yr -1 errors. 



4. The influence of radial velocity 

4.1. Principle of the method 

To quantify the influence of radial velocity on the HTPM proper 
motion for a star, we take the Hipparcos astrometric data (with 
epoch 1991.25) and the literature radial velocity from the XHIP 
catalogue (Section[2]- see Section |o"2"|f or stars without literature 
radial velocity) and use Equations (|6J-((7]l to predict the star's 
celestial position in 2014.5r| i.e., the mean epoch of the first 
intermediate-release Gaia astrometry. We then use backward so- 
lution, i.e., apply Equations ([8]m9]), to recover the proper mo- 
tion from the given Hipparcos and Gaia positions on the sky, as- 
suming that the parallax and time interval are known. Clearly, if 
the radial velocity (radial proper motion) used in the backward 
solution is identical to the radial velocity used in the forward 
propagation, the derived (HTPM) proper motion is essentially 
identical to the input (XHIP) proper motion (Table [TJ. However, 
by varying the radial velocity used for the backward solution 
away from the input value, the sensitivity of the HTPM proper 
motion on radial velocity is readily established. This sensitiv- 
ity does not depend on the sign but only on the magnitude of 



6 cos 2 6o 



-(Aa*) 3 



(8) 



4 To account for flexibility in the launch schedule of Gaia and to be on 
the safe side, we actually used a time interval of 25 years, i.e., 2016.25. 
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Table 1. Proper-motion errors accumulated over t — 25 yr as 
function of declination due to truncation of the full model up 
to and including first, second, and third-order terms in time for 
an 'extreme', i.e., nearby, fast-moving star: m = 500 mas (r 
2 pc), n a — Us — 2000 mas yr , and v r = 50 km s _1 . The unit 
nas stands for nano-arcsecond. 
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Fig. 2. Schematic diagram showing how to quantify the sensitiv- 
ity of the proper motion to a change in (read: measurement error 
of) radial velocity. A change in the radial velocity Av r introduced 
before the backward solution leads to a certain (HTPM) proper- 
motion error Aft. The linear dependence is commented on in 
Section 4.3 Since the magnitude of the proper-motion error does 



not depend on the sign but only on the magnitude of the radial- 
velocity variation, the sensitivity curve is symmetric with respect 
to the true radial velocity. The dashed horizontal line denotes the 
maximum perspective-acceleration-induced proper-motion error 
we are willing to accept in the HTPM proper motion. The dis- 
tance E between the intersection points of the dashed horizontal 
line and the solid sensitivity curves determines the tolerance on 
the radial-velocity error. 




0.4 -0.2 0.2 0.4 
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Fig. 3. Sensitivity of the HTPM proper motion in right ascen- 
sion to radial velocity for HIP70890 (Proxima Centauri). The 
sensitivity is linear a nd h as a value C a = -74.59 jt/as yr~' 
per km s _1 (Section 4.3 1. The dashed horizontal line indi- 



cates the maximum-tolerable perspective-acceleration-induced 
proper-motion error caused by an incorrect radial velocity. 
Since the expected HTPM standard error in right ascension is 
97 yuas yr 1 for this star, we set this threshold to 97/10 = 
9.7 fias yr -1 . This implies, for a confidence level c = 68.27%, 
that the maximum-acceptable radial-velocity standard error cr Vr 
for this object is | • E = 0.13 km s _1 . 



the HTPM proper motion derived from the backward solution. 
The dashed horizontal line represents the maximum perspective- 
acceleration-induced proper-motion error that we are willing to 
accept. If we denote the radial-velocity interval spanned by the 
intersections between the dashed, horizontal threshold line and 
the two solid sensitivity curves by E (either E ff or E 5 ), the toler- 
ance on the radial-velocity standard error is easily expressed as 
-i • E < Av, < j • E. The question now is: what is the prob- 
ability that the error in radial velocity (i.e., true radial velocity 
minus catalogue value) is smaller than - 1 • E or larger than 5 • E? 
Naturally, we want this probability to be smaller than a chosen 
threshold 1 -c, where c denotes the confidence level (for instance 
c = 0.6827 for a' la result'): 



2 A 1 

■ - < Av, < - 
2 2 



2 • CTy 



> c, 



(10) 



where we have assumed that the error distribution for AV r is a 
normal distribution with zero mean and variance crj and where 

<b(x) denotes the error function with argument x / V2: 



O(jc) 



V2 r 
Jo 



exp 



df = erf 



V2 

V^r Jo 

From Equation (fTO), one can easily deduce: 



2 



(11) 



the radial-velocity variation. Figure [2] schematically shows this 
idea. The abscissa shows the change in radial velocity, Av r , with 
respect to the input value used for the forward propagation. The 
ordinate shows Afi (either Afj. a or A/j ), i.e., the difference be- 
tween the input value of the proper motion (either fi a or fig) and 



0-v, < 



2-<D- 1 (c) 



(12) 



where <J> 1 denotes the inverse of O (e.g., O 1 [0.6827] = 1). For 
the 'special case' c = 68.27%, Equation ( 12 1 hence simplifies to 

(T Vr < \ 



1 ■ E. 
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Fig. 4. Histogram of the distribution of £ in the N = 10,000 
Monte-Carlo simulations for star HIP38 (vr = 23.64 ± 0.66 mas, 
so 3% relative error). The smooth, red curve is a Gaussian fit of 
the histogram; it provides a good representation. 
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Fig. 5. Histogram of the distribution of 2 in the N = 10,000 
Monte-Carlo simulations for star HIP8 (m = 4.98 + 1 .85 mas, 
so 37% relative error). The smooth, red curve is a Gaussian fit 
of the histogram; it provides a poor representation and does not 
account for the tail in the distribution. 



4.2. Example for a real star 

Figure [3] is similar to the schematic Figure [2] but shows data 
for a real star, HIP70890. This object, also known as Proxima 
Centauri (or a Cen C), is a nearby (vr = 772.33 ± 2.60 mas), 
fast-moving iji a = -3775.64 + 1.52 mas yr _1 and fig = 768.16 + 
1.82 mas yr~') M6Ve flaring emission-line star which is known 
to have a significant perspective acceleration. HIP70890 is actu- 
ally one of the 21 stars in the Hipparcos catalogue for which the 



perspective acceleration was taken into account ( |ESA 1997a 
Section 1.2.8). Figure [3] was constructed by comparing the in- 
put proper motion used in the forward propagation with the 
proper motion resulting from the backward solution while using 
a progressively differing radial velocity in the backward solution 
from the (fixed) value used in the forward propagation. The ac- 
tual radial-velocity variation probed in this figure is small, only 
±0.5 km s 1 . 

HIP70890 is relatively faint {Hp = 10.7613 mag) and the 
expected HTPM proper-motion standard error is 97 fias yr _I 
(see Section 6.1 for details). If using a ten times lower thresh- 
old, i.e., 9.7 yuas yr~', for the perspective-acceleration-induced 
proper- moti on error caused by an incorrect radial velocity (see 
Section 6.1 for details), we find that 2 = 0.26 km s _1 . This im- 
plies, for a confidence level c = 68.27%, that the radial velocity 
should have been measured for this object with a standard error 
smaller than <x., < - ~ , £ „ = ~ • E = 0.13 km s _1 . The lit- 

2-U)- 1 (r =0.6827) 2 

erature radial velocity for this star is v r = -22.40 + 0.50 km s _I 
(with quality grade 'B'; Section 2.2 1, which is not precise 



enough. New spectroscopic measurements are thus needed for 
this object to reduce the standard error by a factor ~4. 

The discussion above has implicitly focused on the right- 
ascension proper-motion component fi a , and the associated 
since the sensitivty of fig is a factor ~4 less stringent for this 
star. It is generally sufficient to consider the most constraining 
case for a given star, i.e., either or "Eg. Therefore, we drop 
from here on the subscript a and 6 on 2, implicitly meaning 
that it either refers to Y, a or "Zg, depending on which one is 
largest. Typically, this is the largest proper-motion component, 
i.e., \fi a \ > -> X = E a and \p a \ < \p. s \ S = S 5 . 



4.3. Derivation of the sensitivity 

Figure [3] shows that the sensitivity of proper motion to radial 
velocity is linear. This can be understood by substitution of 
Equations |8]l-((9]l in Equations (RjJ-lfTj, after replacing v r , as 
used in the forward propagation, by v r + Av r in the backward 
solution: 



A/j. a = 



Act f Aa |/i r + — )/ - ian<>i,A<i A<)+ 



3 cos 2 £ - 1 
6 cos 2 6q 
1 



(Aa*Y - tan SoAa*A6 1 /J. r + — jf 



a * Av r 

H a t + Aa 1 

r 



Aa* 
r 



tan5oAff*A5 • 



1 



Av, 

• tan^oAff A5 1 

r 



■Av r , 



(13) 



c a 



which immediately shows Afi a oc Av r . A similar analysis for the 
sensitivity coefficient Cg yields: 



Afi s 



A6 1 . , 1 

— + -tan5 (Aa f ■ - 
r 1 r 



■Av r . 



(14) 



The coefficients C a and Cg quantify the proper-motion- 
estimation error caused by a biased knowledge of v r given mea- 
sured displacements Aa* and A6. They can hence more formally 
be defined as the partial derivatives of fi a and /ug from Equations 
<|8j and |9| with respect to v r with Aa* and A6 kept constant: 



C a — 



Cg = 



d\\ 
dpg 
dv r 



dpg 

Vr 



[Aa* - tan<5oAa*A(5] 

1 . 9 

A6 + - tan 5 (Aa ) 



(15) 
(16) 



Clearly, our relations confirm the well-known, classical re- 
sult (e.g., Dravins et al. 1999 Equation 4) that perspective- 
acceleration-induced proper-motion errors are proportional to 
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the product of the time interval, the parallax m oc r _1 , the 
proper-motion components /j. a g, and the radial-velocity error 
Av r : Afi a ,s = C a j ■ Av r = ■ ix a jt ■ Av r . Equations (l3]i-([l4"} 
show that the perspective-acceleration-induced proper-motion 
error caused by a radial-velocity error does not depend on the 
radial velocity v r itself but only on the error Av r . This may look 
counter-intuitive at first sight, since the proper motion itself is 
sensitive to the precise value of the radial velocity. The error 
in the proper motion, however, is sensitive only to the radial- 
velocity error. In other words, the slopes of the V-shaped wedge 
in Figure[3]do not depend on the absolute but only on the relative 
labeling of the abscissa. 



4.4. Taking measurement errors into account 

So far, we ignored the measurement errors of the astrometric pa- 
rameters a, 6, m, jj. a , and A natural way to take these errors 
into account is by Monte-Carlo simulations: rather than deriving 
E once, namely based on the astrometric parameters contained in 
the XHIP catalogue, we calculate E a large number of times (typ- 
ically N = 10,000), where in each run we do not use the cata- 
logue astrometry but randomly distorted values drawn from nor- 
mal distributions centred on the measured astrometry and with 
standard deviations equal to the standard errors of the astromet- 
ric parameters (denoted A^(mean, variance)). We also randomly 
draw the radial velocity in each run from the normal probability 
distribution N(v r , cr 2 ). 

The Monte-Carlo simulations yield N = 10,000 values for 
E; the interpretation of this distribution will be addressed in 
Section [5] Two representative examples of the distribution of E 
are shown in Figures [4] and [5] The first distribution (Figure [4] 
representing HIP38) is for a nearby star with a well-determined 
parallax: m = 23.64 + 0.66 mas, i.e., 3% relative parallax error. 
The smooth, red curve in Figure [4] is a Gaussian fit of the his- 
togram; it provides a good representation of the data. The sec- 
ond distribution (Figure [5] representing HIP8) is for a distant 
star with a less well-determined parallax: m = 4.98 ± 1.85 mas, 
i.e., 37% relative parallax error. This results in an asymmetric 
distribution of E with a tail towards large E values. This is eas- 
ily explained since we essentially have E oc Av, oc C~K oc m~ l , 
meaning that the distribution of E reflects the probability dis- 
tribution function of m~ 



r. The latter is well-known (e.g. 



Kovalevsky 1998; Arenou & Luri| [1999[ ) for its extended tail 
towards large distances and its (light) contraction for small dis- 
tances. The smooth, red curve in Figure[5]is a Gaussian fit of the 
histogram; it provides an inadequate representation of the data. 
We will come back to this in Section IBT21 

To avoid dealing with (a significant number of) negative 
parallaxes in the Monte-Carlo simulations, we decided to ig- 
nore 11,171 entries with insignificant parallax measurements in 
the XHIP catalogue; these include 3,920 entries with m < 
and 7,251 entries with < Tnjcr m < 1 (recall that negative 
parallaxes are a natural outcome of the Hipparcos astromet- 
ric data reduction, e.g., Arenou et al. 1995 i. This choice does 
not influence the main conclusions of this paper: perspective- 
acceleration-induced HTPM proper-motion errors are significant 
only for nearby stars whereas negative and low-significance par- 
allax measurements generally indicate large distances. 



5. Evaluation criteria 

From the Monte-Carlo distribution of E (Section |4~4] i, we want to 
extract information to decide whether the radial velocity avail- 



able in the literature is sufficiently precise or not. For this, we 
develop two evaluation criteria. 

5.1. The Gaussian evaluation criterion 

The first criterion, which we refer to as the Gaussian criterion, 
is based on Gaussian interpretations of probability distributions. 
It can be applied to all stars but, since we have seen that distant 
stars do not have a Gaussian E distribution, but rather a distri- 
bution with tails towards large E values (Figure [5]), the Gaussian 
criterion is unbiased, and hence useful, only for nearby stars. For 
these stars, the Monte-Carlo distribution of E is well described 
by a Gaussian function with mean ji/y. and standard deviation cry, 
(Figure [4|. The equations derived in Section 4.1 by ignoring as- 
trometric errors are easily generalised by recognising that both 
the radial-velocity error and the distribution of E have Gaussian 
distributions, with standard deviations cr,, and crj_, respectively: 



P|-| <Av,.< || = d) 



2 • <Ty t 



> c 



generalises to: 



/ 



dE 



4 



2ncr\ 



exp 



2 • <r v 



> c 



(17) 



(18) 



Using Equation 8.259.1 from Gradshteyn & Ryzhik ( 2007) , 
Equation ( 18 i simplifies to: 



Ms 



4 



4- cr? +oi 



> c, 



(19) 



which correctly reduces to Equation ( fT0| ) for the limiting, 'error- 
free' case yUj; — > E and crj_ — > in which the Gaussian distribu- 
tion of E collapses into a delta function at y^. = E, i.e., 5(1,). 

For the example of HIP70890 discussed in Section |4.2| we 
find ^ 2 = 0.26 km s -1 and cr y = 0.000, 89 km s _I so that, with 
cr Vr = 0.50 km s" 1 , Equation |l9| returns c = 20.58%. 



It is trivial, after re-arranging Equation ( 19 1 to: 



ov < 



m! 



crt 



V4-[0-!(c)] 2 4 



(20) 



to compute the required standard error of the radial veloc- 
ity to comply with a certain confidence level c. For instance, 
if we require a c = 99.73% confidence level (for a '3cr re- 
sult'), the radial-velocity standard error of HIP70890 has to be 
0.04 km s _1 . We finally note that Equation ( 20 1 correctly reduces 
to Equation ( 12 1 for the limiting, 'error- free' case — > E and 
cr 2 — > 0. 



A limitation of Equation ( 20 » is that the argument of the 
square root has to be non-negative. This is physically easy to 
understand when realising that, in the Gaussian approximation, 
one has cr-^ ss ■ (cr w jTn) (see Section 4.4), so that: 



Ml 



> -> _1 (c) < 



4-[0- 1 (c)] 2 4 



c < $ 



■ (21) 



So, for instance, if a certain star has cr/ov = 2 (a '2cr parallax'), 
the Gaussian methodology will only allow to derive the radial- 
velocity standard error <x, v required to meet a confidence level 
c = 95.45% or lower. 
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Fig. 6. Predicted HTPM proper-motion error as function 
of the Hp Hipparcos broad-band magnitude following 
Equations (23i-(24i. We require perspective-acceleration- 



induced proper-motion errors to be an order of magnitude 
smaller (factor of safety = FoS = 10; see Section 7.2 1. 



5.2. The robust criterion 

Since the Monte-Carlo distribution of S values is not Gaussianly 
distributed for distant stars (Figure [5]), the Gaussian criterion re- 
turns incorrect estimates; in fact, the estimates are not just incor- 
rect but also biased since the Gaussian criterion systematically 
underestimates the mean value of 2 (i.e., p£) and hence sys- 
tematically provides too conservative (small) estimates for cr Vr 
through Equation (20 1. Rather than fitting a Gaussian function, 



we need a more robust estimator of the location and width of the 
2 distribution than the Gaussian mean and standard deviation. 
This estimator is contained in the data itself and provides, what 
we call, the robust criterion. 

Let us denote the individual values of 2 derived from the 
N = 10,000 Monte-Carlo simulations by 2,, with i = l,...,N. 
Equation ( 18 1 is then readily generalised for arbitrary distribu- 
tions of 2 to: 



— < Av r < — 

2 2 



\ N 1 

) ^ 



2 ■ cr„ 



> c. 



(22) 



The inverse relation generalising Equation ( |20[ > by expressing 
cr,, as function of c, required to determine the precision of the 
radial velocity required to comply with a certain confidence level 
c, is not analytical; we hence solve for it numerically. 

The robust criterion generalises the Gaussian criterion. Both 
criteria return the same results for nearby stars which have a 
symmetric Gaussian distribution of S values. In general, there- 
fore, the robust criterion is the prefered criterion for all stars, 
regardless of their distance. 



6. Application of the evaluation criteria 

6.1. Target proper-motion-error threshold 

Before we can apply the Gaussian and robust criteria, we 
have to decide on a target proper-motion-error threshold for 
each star (i.e., the location of the dashed horizontal lines in 
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Fig. 7. The fraction of stars with XHIP literature radial velocities 
which are contained in the radial-velocity interval [m - R, m + R] 
as function of R, with m - -2.21 km s the mean v r for the 
Gaussian criterion and m = -2.00 km s _1 the median v, for the 
robust criterion (Sections 12.21 and 16. 2b. For the Gaussian crite- 



rion, we represent the histogram of literature radial velocities by 
a Gauss with standard deviation cr = 22.44 km s _1 (Figure nil. 
The dashed lines represent the classical limits Icr = 68.27%, 
2cr = 95.45%, and 3cr = 99.73%. The fraction of stars with 
the robust criterion builds up more slowly as a result of the non- 
Gaussian broad wings as well as outliers representing halo and 
runaway stars. Since the Gaussian criterion ignores these fea- 
tures, it returns biased results for stars without literature radial 



velocity (see Section 6.3 i 



Figures [2] and 13}. We adopt as a general rule that the max- 
imum perspective-acceleration-induced HTPM proper-motion 
error caused by radial-velocity errors shall be an order of mag- 
nitude smaller than the predicted standard error of the HTPM 
proper motion itself (Section 7.2 discusses this choice in more 
detail). The latter quantity has been studied by Mignard (2009) 
and can be parametrised as: 



cr Ma [//as yr 1 ] = -227.8 + 122.1 • H - 20.39 • H 2 

+ 1.407 H 3 -0.02841 ■ H 4 
cr Ms [juasyr -1 ] = 127.2 - 47.0 H + 8.30 ■ H 2 



-0.686 H 3 +0.02581 ■ H 4 , 



(23) 
(24) 



where H = max{6.5, Hp [mag]} with Hp the Hipparcos broad- 
band magnitude. These relations are shown in Figure|6] The pre- 
dicted HTPM standard errors include residual errors caused by 
the correction for the parallactic effect in the Gaia data (see foot- 
note 1), the expected number and temporal distribution of the 
Gaia field-of-view transits for the Gaia nominal sky scanning 
law, and the expected location-estimation precision fcentroid- 
ing error') of Gaia's CCD-level data. 

6.2. Stars without literature radial velocities 



As already discussed in Section 
XHIP catalogue do not have a 



2.2[ the majority of stars in the 
iterature radial velocity. These 



7 1 , 563 objects are treated as stars with radial velocity, with three 
exceptions: 
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Fig. 8. Comparison of the robust and Gaussian criteria for all 
stars with cr w l%n better than 5%. The top panel compares the 
confidence levels while the bottom panel compares the radial- 
velocity standard errors cr,, required to reach a confidence level 
c = 68.27%. The top panel shows two branches of data points: 
the linear, one-to-one branch corresponds to stars with a mea- 
sured radial velocity, whereas the lower, curved branch corre- 
sponds to stars without measured radial velocity in the XHIP 
catalogue. As explained in Section 6.2 the latter objects suffer 
from a bias in the Gaussian confidence level cc aU ss- 



1. For each of the N = 10,000 Monte-Carlo simulations, the 
radial velocity is randomly taken from the full list of 46, 392 
radial velocites contained in the XHIP catalogue (Figure [T}. 
In practice, this choice does not influence our results since 
we are not sensitive to the absolute value of the radial veloc- 
ity (Section |4.3| >. But at least in principle this choice means 
that there is a finite probability to assign a halo-star-like, i.e., 
large radial velocity in one (or more) of the Monte-Carlo 
runs. Regarding the HTPM Project, the best choice for stars 
without known radial velocity is to use v r = -2.00 km s _1 , 
which is the median value of the distribution (recall that the 
mean equals -2.21 km s _1 ). 

2. For the Gaussian evaluation criterion (Section |5.1| l, we use 
<x, v = 22.44 km s _1 from the Gaussian fit to the distribution 
of all literature radial velocities. We also follow this recipe 
for the Gaussian criterion for the 1,753 entries which do 
have a radial velocity but which do not have an associated 
standard error in the XHIP catalogue (this concerns 23 en- 
tries with quality grade 'C and 1,730 entries with quality 
grade 'D'). Clearly, this approach ignores the broad wings of 
the distribution visible in Figure [T] as well as a small but fi- 
nite number of halo stars and runaway stars with heliocentric 
radial velocities up to plus-or-minus several hundred km s" 1 
(see Section |2~2] >. As a result, the Gaussian criterion system- 
atically returns an overly optimistic (i.e., too large) confi- 
dence level CGauss for stars without literature radial velocity 
(Figures [7] and [8]). This bias comes in addition to the bias 



for distant stars for which the Gaussian criterion returns too 
conservative (small) estimates for cr,, (Section 5.2 1. 
For the robust evaluation criterion (Section |5.2) , we do not 
make a priori assumptions except that the overall radial- 
velocity distribution of stars without known radial velocity 
is the same as the distribution of stars with literature radial 
velocity, including broad wings and 'outliers'. We thus use 
Equation ( 22 1 in the form: 



Z 



1 / Si Si 
-■PI — - < Av r < — 
N \ 2 2 



> c, 



(25) 



where the probability P that Av, is contained in the interval 
[—5 ■ Sj, I ■ Sj] is calculated as the fraction of all stars with 
literature radial velocities in the XHIP catalogue which has 
v r e [-| • Sj + v r , me dian, 3 • Sj + vidian] (recall that the me- 
dian radial velocity equals -2.00 km s _1 ). We thus cater for 
the broad wings of the observed radial-velocity distribution 
(Figures[TJand|7| as well as for the probability that the object 
is a (fast-moving) halo star, avoiding the bias in the Gaussian 
criterion discussed in the previous bullet. 



6.3. Results of the application 

We applied the Gaussian and robust criteria, as described in 
Sections 5.1 and 5.2 respectively, to the XHIP catalogue. The 
results are presented in Table |A. 1~| (Appendix [A]), which is avail- 
able electronically only. The run time for N = 10,000 Monte- 
Carlo simulations is typically ~0.7 s per star and processing the 
full set of 117,955 - 11,171 = 106, 784 XHIP entries with sig- 
nificant parallaxes (Section |4~4| hence takes about one day. 

Figures |9"HT0] and Tables |2J-[3] show the results for the confi- 
dence levels of the literature radial velocities contained in the 
XHIP catalogue. We find, not surprisingly since perspective- 
acceleration-induced proper-motion errors are relevant only for 
nearby, fast-moving stars - which are relatively rare - that the 
majority of stars have confidence levels exceeding c = 99.73%. 
This indicates that, at the c = 99.73% confidence level, the avail- 
able radial velocity is sufficiently precise or, for stars without 
literature radial velocity, that the absence of a literature radial 
velocity, and hence the assumption v, = -2.00 km s _1 for the 
robust criterion or v r = -2.21 km s for the Gaussian criterion, 
is acceptable. This holds for more than 100,000 stars using the 
robust criterion (Table [3j and more than 85,000 stars using the 
Gaussian criterion (Table |2j. The large difference between the 
two criteria does not come unexpectedly: 



1. 



We already argued in Section [572] that the Gaussian criterion 
is biased for distant stars, with distant meaning that the par- 
allax probability distribution has an associated asymmetric 
distance probability distribution. For such stars, the Gaussian 
criterion systematically underestimates the mean value of S 
and hence returns too conservative (small) values for o~ Vr for 
a given value of the confidence level c and too pessimistic 
(small) values of c for a given value of <x,, r ; 
We already argued in Section 6.2 that the Gaussian criterion 



is biased for stars without literature radial velocity. For such 
stars, the Gaussian criterion systematically returns too opti- 
mistic (large) values of the confidence level c since it ignores 
the broad wings of the observed distribution of radial veloci- 
ties (Figures [T] [7] and[8]l and also ignores the probability that 
the object is actually a halo (or runaway) star. 
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Table 2. Number of stars as function of confidence level c es- 



tablished using the Gaussian criterion from Section 5.1 The 
number in brackets in column 2 indicates the number of bright 
stars (G < 5.7 mag or Hp < 6.0 mag), i.e., those stars not de- 
tectable with Gaia and hence not contained in the HTPM cat- 
alogue (Section [7. 1| >. Stars with insignificant parallax measure- 
ments have not been processed and the total number of entries 
hence equals 117,955 - 11, 171 = 106,784 (see Section[44i. 



Conf. level c [%] 


# stars (# bright) 


v, e XHIP 


v r <t XHIP 


[0 - 68.27) 


225 ( 12) 


108 


117 


[68.27 - 95.45) 


2,087 ( 27) 


598 


1,489 


[95.45 - 99.73) 


18,815 ( 88) 


4,104 


14,711 


[99.73 - 100] 


85,657 (4, 341) 


38,700 


46, 957 


Total 


106,784(4, 468) 


43,510 


63,274 



Table 3. Number of stars as function of confidence level c es 

The num 



tablished using the robust criterion from Section 5.2 
ber in brackets in column 2 indicates the number of bright stars 
(G < 5.7 mag or Hp < 6.0 mag), i.e., those stars not de- 
tectable with Gaia and hence not contained in the HTPM cat- 
alogue (Section [7. l| i. Stars with insignificant parallax measure- 
ments have not been processed and the total number of entries 
hence equals 117,955 



11,171 = 106,784 (see Section 4.4 1 



Conf. level c [%] 


# stars (# bright) 


v, e XHIP 


v r i XHIP 


[0 - 68.27) 


206 ( 12) 


97 


109 


[68.27 - 95.45) 


1, 112 ( 15) 


150 


962 


[95.45 - 99.73) 


5,244( 10) 


135 


5, 109 


[99.73 - 100] 


100, 222(4,431) 


43, 128 


57,094 


Total 


106,784(4, 468) 


43,510 


63,274 



The robust criterion does not suffer from these biases and hence, 
being more reliable, is prefered for all stars. The Gaussian crite- 
rion, nonetheless, provides a useful and also easily interpretable 
reference test case and we hence decided to retain it. Figure [8] 
shows that, for nearby stars with literature radial velocities, the 
Gaussian and robust criteria return equivalent results. 

For a small but non-negligible number of stars, Table [3] in- 
dicates insatisfactory results: 206 stars have a confidence level 
c < 68.27%: 97 of these do have a literature radial velocity 
in the XHIP catalogue but one which is insufficiently precise. 
The remaining 109 stars do not have a spectroscopically mea- 
sured radial velocity (at least not one contained in the XHIP cat- 
alogue). New spectroscopy is hence required for these stars to 
guarantee a confidence level of at least 68.27%. For increased 
confidence levels, the numbers obviously increase: if requiring 
a c — 99.73% confidence level for all objects, for instance, the 
number of 'problem stars' increases to 6, 562, split into 382 with 
insufficiently-precise known radial velocity and 6, 180 without 
known radial velocity. We conclude that, depending on the confi- 
dence level one wants to achieve, hundreds to thousands of stars 
need to be spectroscopically re-measured. 



Figure 1 1 shows the robust confidence level versus Hp mag- 
nitude. One can see that the typical star which needs a high- 
priority spectroscopic measurement (i.e., c < 68.27%) has Hp in 



the range 8-12 mag. Figure 12 shows the radial-velocity preci- 
sion required to reach c - 68.27% (computed with the robust cri- 
terion) versus magnitude. Precisions vary drastically, from very 




50 100 



100 
Gaussian 



50 100 
confidence level 



50 100 
[%] 



Fig. 9. Histograms of the Gaussian confidence level c from 
Section [B~T| for all stars combined, stars without measured ra- 
dial velocity, and stars with measured radial velocity as func- 



tion of radial-velocity quality grade Q Vi (see Section 2.2 and 
AppendixfA). The vast majority of objects have c > 95.45% (see 
Table |2j; they have been omitted from the histograms to improve 
their legibility. 




100 




100 



50 100 50 100 50 
Robust confidence level [%] 

Fig. 10. As Figure]^ but for the robust criterion from Section 



100 



5.2 



stringent values well below 1 km s 1 to very loose values, up to 
several tens of km s . 



7. Discussion 

7.1. HTPM bright limit 

The HTPM catalogue will contain the intersection of the 
Hipparcos and Gaia catalogues. Whereas the Hipparcos cata- 
logue ( |ESA| [T997b| l, which contains 117,955 entries with as- 
trometry (and 118,218 entries in total), is complete to at least 
V = 7.3 mag ( |ESA| 1 1997a] Section 1.1), the Gaia catalogue 
will be incomplete at the bright end. Gaia's bright limit is G = 
5.7 mag (de Bruijne 2012 1, where G is the white-light, broad- 
band Gaia magnitude, which is linked to the Hipparcos Hp, the 



Cousins /, and the Johnson V magnitudes through (Jordi et al. 

|20T0l ): 



G-V = -0.0447 - 0.1634 • (Hp - I) + 

+0.0331 • (Hp - I) 2 - 0.0371 • (Hp - 1) 3 . 



(26) 
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Hipparcos broad-band Hp magnitude [mag] Hipparcos broad-band Hp magnitude [mag] 



Fig. 11. Robust confidence level c versus Hipparcos broad-band 
Hp magnitude. The red box in the bottom-right comer denotes 
the approximate area for high-priority follow-up spectroscopy: 
stars with confidence level c < 68.27% and Hp > 6 mag. The 
latter restriction roughly reflects Gaia's - and hence HTPM's - 
bright limit (Section [7~T) . 



For stars in the Hipparcos catalogue, the colour G - Hp ranges 
between -0.5 and 0.0 mag, with a mean value of -0.3 mag. This 
means that G = 5.7 mag corresponds roughly to Hp * 6.0 mag. 
In practice, therefore, we are not concerned with the bright- 
est ~4, 509 stars in the sky. The expected number of HTPM- 
catalogue entries is therefore -117, 955-4, 509 ~ 1 13, 500. The 
number of entries with significant parallax measurements equals 
117,955 - 11,171 = 106,784, of which 106,784 - 4,468 = 
102, 316 have G > 5.7 mag. 



7.2. Proper-motion-error threshold 

For both the Gaussian and the robust criteria, we adopt, 
somewhat arbitrarily, the rule that the perspective-acceleration- 
induced HTPM proper-motion error caused by radial-velocity 
errors shall be an order of magnitude smaller than the predicted 



standard error of the HTPM proper motion itself (Section 6. 1 
The adopted Factor of Safety (FoS) is hence 10. Some readers 
may find that this 'rule' is too stringent. Unfortunately, there is 
no easy (linear) way to scale our results if the reader wants to 
adopt a different value for the FoS. Clearly, the value of Z in 
each Monte-Carlo run is linearly proportional to the FoS (see, 
for instance, Figure [2]). However, the robust confidence level 
Crobust and the radial-velocity standard error cr Vr required to meet 
a certain value of c ro b us t do not linearly depend on Z but on 
the properties of the, in general asymmetric, distribution of the 
N = 10,000 values of Z resulting from the Monte-Carlo pro- 
cessing. To zeroth order, however, one can assume a linear re- 
lationship between cr Vr and Z and hence the adopted FoS, as is 
also apparent from the 'error-free' criterion cr,, < | • X derived 
for c = 68.27% in Section [4~T] This is in particular a fair approx- 
imation for small variations around the default value (FoS = 10) 
in combination with nearby stars, which are most interesting be- 
cause these are most sensitive to perspective acceleration. The 
Monte-Carlo distribution of Z values for these objects is gener- 
ally well behaved, i.e., symmetric and with cr m <k m and hence 
o"z <k py. (see, for instance, Figure |4j see also Section |5. \\ . 



Fig. 12. Radial-velocity precision (standard error) for stars with 
robust confidence level c < 68.27% required to upgrade their 
confidence level to c — 68.27%. 



Figure [13] shows, as an example for the nearby star HIP57367 
(see Table]?]), how <r Vr (required for c 10 b U st = 68.27%) and c r0 b us t 
vary as function of the FoS. Linear scaling around the default 
FoS = 10 provides a decent approximation, at least over the 
range 10/2 = 5 < FoS < 20 = 10 • 2. 



7.3. Number of Monte-Carlo simulations 

To take measurement errors in the XHIP astrometry and liter- 
ature radial velocities into account, we adopt a Monte-Carlo 
scheme in which we run a number N of Monte-Carlo simula- 
tions for each star in which we randomly vary the astrometric 
and spectroscopic data within their respective error bars (see 
Section 4.4 1. Clearly, the higher the value of N is, the more re- 
liable the results are. We adopted N = 10,000 as a practical 
compromise, resulting in an acceptable, typical run time of ~1 s 
per star as well as smooth distributions of Z (see, for instance, 
Figure [4] or |5j. To investigate the repeatability and hence relia- 
bility of our robust confidence levels c ro b us t and radial-velocity 
errors o~ Vt , we have repeated the entire processing with N = 
10,000 runs 100 times for the 206 stars with c rob ust < 68.27% 
(Section 6.3 1 and find that the typical variation of the confidence 



level and the radial-velocity error, quantified by the standard de- 
viation divided by the average of the distribution containing the 
100 results, is less than 0.1% and 0.2%, respectively; the max- 
imum variation among the 206 objects is found for HIP107711 
and amounts to 0.5% and 0.6%, respectively. 



7.4. High-priority and challenging stars 

Table |4] shows the ten stars with the lowest robust confidence 
level. These stars are the highest-priority targets for spectro- 
scopic follow-up. Nine of the ten entries do have a literature ra- 
dial velocity but one which is insufficiently precise. The highest- 
priority object (HIP57367) does not yet have a spectroscopic 
radial velocity and needs a measurement with a standard error 
better than 1 km s^ 1 . For this particular object, this challenge 
seems insurmountable since it is one of the 20 white dwarfs with 
Hipparcos astrometry (Vaucl air et al.| [l997), objects for which it 
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5 10 15 

Factor of Safety [FoS] 



Fig. 13. Dependence, for HIP57367 (see Table |4|, of cr Vr (re- 
quired to reach c ro b ust = 68.27%) and the robust confidence 
level c ro bust on the Factor of Safety (FoS), i.e., the mini- 
mum factor between the predicted HTPM standard error and 
the perspective-acceleration-induced HTPM proper-motion er- 
ror caused by radial-velocity errors. The default FoS value 
adopted in this study is 10. The object is representative for a 
nearby star with a well-determined parallax. The straight, red 
lines indicate linear scaling relations starting from the default 
value FoS = 10. 



7.5. Object-by-object analyses and other literature sources 

In general, and in particular for the most interesting, delicate, or 
border cases, it will be useful to perform a more in-depth litera- 
ture search for and study of radial velocities and other available 
data before embarking on ground-based spectroscopy. For in- 
stance, we found a SIMBAD note on the Hipparcos catalogue 
( |ESAl|1997a|bl CDS catalogue 1/239) for HIP114110 (c robust = 
71.23%) and HIP1 14176 (c rob ust = 60.00%) that they are non- 
existing objects: "HIPl 141 10 (observed with ///P14113^} and 
HIPl 14176 (observed with ///PI 14177) are noted as proba- 
ble measurements of scattered light from a nearby bright star. 
The non-reality of 114110 and 114176 (traced to fictitious en- 
tries in the WDS and INC A ) has been confirmed by MMT ob- 
servations reported by D. Latham (private communication, 8 
May 1998), and confirmed by inspection of the DSS [J.L. Falin, 
12 May 1998]". In addition, the completeness and coverage 
level of the XHIP radial-velocity compilation is not known. We 
did query SIMBAD as well as the Geneva-Copenhagen-Survey 
(GCS, CORAVEL) database flNordstrom etal] [20041 CDS cata- 
logue V/l 17) and the RAVE database ( |Siebert et al.[|20TT] CDS 
catalogue III/265, with 77, 461 entries with a mean precision of 
2.3 km s ) for radial velocities for the 109 stars without XHIP 
radial velocity and with confidence level below 68.27% but did 
not find new data. Unfortunately, the treasure contained in the 
full CORAVEL database (45, 263 late-type Hipparcos stars with 
precisions below 1 km s _1 ), the public release of which was 
announced in |Udry et al. ( 1997 1 to be before the turn of the pre- 
vious millennium, remains a mystery to date. All in all, dedi- 
cated studies for individual objects might pay off by reducing 
the needs for spectroscopic follow-up. 



is notoriously difficult to obtain - even low-precision - spectro- 
scopic radial velocities. 

Table [5] shows the ten stars, among the subset of stars with 
unacceptably-low robust confidence level (c ro b us t < 68.27%), 
with the smallest radial-velocity standard errors required to raise 
the robust confidence level to c ro b us t = 68.27%. Since c ro b us t < 
68.27%, these stars do clearly need spectroscopic follow-up. 
However, the radial-velocity standard errors reach values as 
small as 0.04 km s , which is a real challenge, not only in terms 
of the required signal-to-noise ratio of the spectroscopic data but 
also in view of the definition of the radial-velocity zero-point at 
this level of precision (Crifo et al. 20101 as well as potential 
systematic errors in the radial velocities, both with instrumental 
origin and with astrophysical causes such as radial- velocity dif- 
ferences between various absorption lines etc. (seelLindegren & 



7.6. Urgency of the spectroscopic follow-up 



Dravins 2003 for a detailed discussion of this and other effects). 



Ta ble 2 in |Dravins et al. (19 99 ) shows the top-39 of stars 
in the Gliese & Jahreiss ( 1995) preliminary third catalogue of 
nearby stars ranked according to the magnitude of the perspec- 
tive acceleration (which is propertional to m ■ p). Similarly, 



Table 1.2.3 in ESA] p997a) shows the top-21 of stars in the 
Hipparcos catalogue for which the magnitude of the perspec- 
tive acceleration is significant enough to have been taken into 
account in the Hipparcos data processing (the accumulated ef- 
fect on position is proportional to w-p ■ |v, |). On the contrary, the 
top- 10 Tables [4] and [5] have been constructed based on the sen- 
sitivity of the perspective acceleration to radial-velocity errors 
(Section pOj ) and the associated confidence level of the available 
literature radial velocity. Hence, although there is a significant 
overlap of stars between the various tables, they are understand- 
ably not identical. 



Mign ard| ( |2009| l already acknowledges that, since the 
perspective-acceleration-induced proper-motion error can 
be calculated as function of radial velocity, a factor - effectively 
the sensitivity coefficients C a and Cg from Equations (fI3}-([T4"|i 



in Section 4.3 



can be published to correct the HTPM proper 
motion for a particular star a posteriori when v, becomes known 
or when a more precise v, becomes available. Therefore, both 
the reference radial velocity v r and parallax m used in the HTPM 
derivation will be published together with the proper-motion 
values themselves. This means that the spectroscopic follow-up 
identified in this paper is not time-critical: the HTPM catalogue 
can and will be published in any case, even if not all required 
spectroscopic follow-up has been completed. Of course, the 
implication for stars without the required radial-velocity knowl- 
edge will be that their HTPM proper motions will include a 
(potentially) significant perspective-acceleration-induced error. 



8. Conclusions 

We have conducted a study of the requirements for the availabil- 
ity of radial velocities for the Hundred-Thousand-Proper-Motion 
(HTPM) project (Mignard, 2009). This unique project will com- 
bine Hipparcos astrometry from 1991.25 with early-release Gaia 
astrometry (~2014.5) to derive long-time-baseline and hence 
precise proper motions. For the nearest, fast-moving stars, the 
perspective acceleration of the objects on the sky requires the 
presence of radial velocities for the derivation of the proper mo- 
tions. We have quantitatively determined, for each star in the 



This is a typo and must be HIPl 141 13. 
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Table 4. The ten stars with the lowest robust confidence level c ro bust m the XHIP catalogue. The column with header 'cr Vt , c ro b us t = 
68.27%' denotes the radial-velocity standard error that should be targeted in the spectroscopic follow-up to raise the confidence 
level of these stars to c ro b us t = 68.27%. The table contains two white dwarfs: HIP57367 (L145 - 141, type DQ, i.e., with carbon 
absorption features) and HIP3829 (van Maanen 2, type DZ, i.e., with metal absorption features). 



HIP 


Vr 


CT Vr 


C robust 


0">y, Crobust = 68.27% 


m 




Hp 


Pa 


Ps 




[km s" 1 ] 


[km s- 1 ] 


[%] 


[kms- 1 ] 


[mas] 


[mas] 


[mag] 


[mas yr~'] 


[mas yr -1 ] 


57367 






2.85 


0.85 


217.01 


2.40 


11.5851 


2664.98 


-348.60 


73182 


35.63 


10.06 


4.22 


0.53 


168.77 


21.54 


8.1824 


961.78 


-1677.83 


86990 


-115.00 


21.00 


4.70 


1.24 


171.48 


2.31 


10.7574 


-1119.35 


-1352.81 


86214 


-60.00 


21.00 


6.16 


1.62 


196.90 


2.15 


10.8964 


-708.98 


-937.40 


3829 


263.00 


4.90 


12.90 


0.80 


234.60 


5.90 


12.5592 


1236.90 


-2709.19 


72511 


-39.60 


10.00 


13.40 


1.65 


235.24 


44.85 


11.8102 


-1389.70 


135.76 


113229 


46.60 


10.00 


16.51 


2.08 


116.07 


1.19 


10.3732 


-1027.76 


-1060.78 


72509 


-38.80 


10.00 


16.87 


2.07 


214.67 


43.88 


12.3140 


-1416.49 


-270.45 


80018 


46.60 


10.00 


18.40 


2.33 


118.03 


2.52 


10.5964 


-740.10 


997.40 


55042 


-35.00 


10.00 


20.15 


2.55 


78.91 


2.60 


11.6071 


-2466.98 


1180.09 



Table 5. The ten stars, among the subset of stars with unacceptably-low robust confidence level (c ro bust < 68.27%), with the most 
stringent radial-velocity-error requirements <x Vf needed to reach c ro b us t = 68.27%. 



HIP 


C~v r ^robust 


0">y, C robust = 68.27% TD CT m 


Hp 


Pa 


Pi 




[kms" 1 ] [kms" 1 ] [%] 


[kms -1 ] [mas] [mas] 


[mag] 


[mas yr -1 ] 


[mas yr _I ] 



87937 


-110.51 


0.10 


27.84 


0.04 


548.31 


1.51 


9.4901 


-798.80 


10277.30 


54035 


-84.69 


0.10 


51.74 


0.07 


392.64 


0.67 


7.5062 


-577.00 


-4761.80 


1 14046 


8.81 


0.10 


59.77 


0.08 


305.26 


0.70 


7.4182 


6768.20 


1327.52 


70890 


-22.40 


0.50 


20.59 


0.13 


772.33 


2.60 


10.7613 


-3775.64 


768.16 


439 


25.38 


0.22 


53.99 


0.16 


230.42 


0.90 


8.6181 


5635.74 


-2338.18 


73182 


35.63 


10.06 


4.22 


0.53 


168.77 


21.54 


8.1824 


961.78 


-1677.83 


84140 


-30.10 


1.10 


52.48 


0.79 


158.17 


5.02 


9.3754 


252.60 


-1571.80 


10138 


57.00 


0.80 


67.40 


0.79 


92.74 


0.32 


6.2571 


2150.30 


673.20 


3829 


263.00 


4.90 


12.90 


0.80 


234.60 


5.90 


12.5592 


1236.90 


-2709.19 


57367 






2.85 


0.85 


217.01 


2.40 


11.5851 


2664.98 


-348.60 



Hipparcos catalogue, the precision of the radial velocity that is 
required to ensure that the perspective-acceleration-induced er- 
ror in the HTPM proper motion caused by the radial-velocity er- 
ror is negligible. Our method takes the Hipparcos measurement 
errors into account and allows the user to specify his/her own 
prefered confidence level (e.g., 68.27%, 95.45%, or 99.73%). 
The results are available in Table A. 1 (Appendix |A} , which 
is published electronically only. We have compared the radial 
velocity-precision requirements to the set of 46, 392 radial ve 



locities contained in the XHIP compilation catalogue (Anderson 



& Francis 2012 1 and find that, depending on the confidence 



level one wants to achieve, hundreds to thousands of stars require 
spectroscopic follow-up. The highest-priority targets are 206 ob- 
jects with a confidence level below 68.27%; 97 of them have a 
known but insufficiently precise radial velocity while the remain- 
ing 109 objects have no literature radial velocity in the XHIP 
compilation catalogue at all. The typical brightness of the objects 
requiring their radial velocity to be (re-)determined is Hp « 8- 
12 mag and the radial-velocity precisions vary drastically, rang- 
ing from 0.04 km s _1 for the most extreme case (HIP87937, also 
known as Barnard's star) to a few tens of km s -1 . With only few 
exceptions, the spectral types are K and M; 73% of them are in 
the south. Gaia's Radial- Velocity Spectrometer (RVS; [Cropper| 
|& Katz| |201 1| > will deliver radial velocities for all stars in the 
HTPM catalogue with Gaia-end-of-mission precisio ns below a 
few km s _1 (and ~10 km s 1 for early- type stars; |de Bruijne| 
|2012[ l; however, these performances require full calibration of 
the instrument and data and hence will most likely only be 



reached in the final Gaia data release, at which time the HTPM 
proper motions will be superseded by the Gaia proper motions. 
Fortunately, the spectroscopic follow-up is not time-critical in 
the sense that the HTPM catalogue will be published with infor- 
mation (sensitivity coefficients and reference parallax and radial 
velocity) to correct the proper motions a posteriori when (im- 
proved) radial velocities become available. We finally note that 
the spectroscopic follow-up requirements for the HTPM proper 
motions quantified in this work will be dwarfed by the require- 
ments coming from the end-of-mission Gaia proper motions, 
to be released around ~2021: for instance for the stars in the 
HTPM catalogue, for which the HTPM proper-motion standard 
errors are 30-190 pas yr 1 , the Gaia proper-motion st andard er- 
rors reach the bright-star floor around 3-4 pas yr" 1 ( |de Bruijne| 
2012), which means that the spectroscopic requirements for the 
correction of perspective acceleration in the Gaia astrometry will 
be a factor ~ 10-50 more demanding. 
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Appendix A: The results data file 



Table [AT] describes the results data file, which is available elec- 
tronically only. 
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Table A.l. Description of the final-results data file. The table, which is only available electronically, covers the 117,955 entries 
of the XHIP catalogue (Section EJ. These are all Hipparcos-catalogue entries with astrometry. Columns 1-7 refer to data extracted 
from the XHIP catalogue described in Section 2] Columns 8-10 summarise some key quantities of our method. Columns 11-14 
provide data from the N = 10,000 Monte-Carlo simulations that we ran for each star (Section [4.4). Columns 15-18 refer to the 



Gaussian criterion developed in Section 5.1 Columns 19-22 refer to the robust criterion, which is described in Section 5.2 



Column 


Value 


Unit 


1 


HIP 


— 


2 


Vr 


km s _1 


3 


< r r, 


km s 


4 


Qvr 




5 


UJ 


mas 


6 




mas 


7 


Hp 


mag 


8 


C 


yuas yr~' per km s 


9 


s 


km s -1 


10 


Flag 


- 


11 


^median 


km s -1 


12 


^smallest 


km s -1 


13 




km s" 1 


14 


""I 


km s -1 


15 


c Gauss 


% 


16 


O" v r ,68.27%,Gauss 


km s' 1 


17 


O" vy,95.45%,Gauss 


km s -1 


18 


C >y.99.73%,Gauss 


km s -1 


19 


'-robust 


% 


20 


C iy,68.27%,robust 


km s -1 


21 


C v r ,95.45%,robust 


km s" 1 


22 


0~ vv,99.73%,robust 


km s -1 



Explanation 



Hipparcos identifier 

Radial velocity in the XHIP catalogue (Section|2.2|l 
Standard error of the radial velocity in the XHIP catalogue (Section 
Quality grade of XHIP radial velocity (1 = A, 2 = B, 3 = C, 4 = 
Parallax from the XHIP catalogue (Section pTl 
Parallax standard error from the XHIP catalogue (Section [2.l[ 
Hipparcos broad-band magnitude 



2.2 1 
D - Section 



2.2 



4.3 i 



Sensitivity of the proper motion to radial velocity (C = d/-i/dv r ; Section 
Value of E for the 'error-free' astrometry from the XHIP catalogue (Section 
Flag indicating whether columns 8-9 refer to right ascension (=1) or declination (= 2) 



4.1 



Median value of the N = 10, 000 Monte-Carlo S values (Section 4.4 

Smallest value of S among the N = 10, 000 Monte-Carlo simulations 

Mean of the Gauss fit to the histogram of the N = 10, 000 Monte-Carlo Z values 

Standard deviation of the Gauss fit to the histogram of the N = 10, 000 Monte-Carlo £ values 



Gaussian confidence level c Gai i SS (Section 5.1 1 

Radial-velocity precision needed for a CGauss = 68.27% confidence level 
Radial- velocity precision needed for a c Gau ss = 95.45% confidence level 
Radial-velocity precision needed for a CQauss = 99.73% confidence level 



5.2i 



Robust confidence level c robllst (Section 
Radial-velocity precision needed for a c, 
Radial-velocity precision needed for a c rob ust 
Radial-velocity precision needed for a c ro t, U st 



: 68.27% confidence level 
: 95.45% confidence level 
: 99.73% confidence level 



Detailed notes: 

- Columns 2-4: for stars without literature radial velocity, these columns list NaN (see Section [2~2) ; 

- Columns 8-22: for the 11, 171 entries with insignificant parallax measurements in the XHIP catalogue, i.e., 3,920 entries with m < and 
7, 251 entries with < nr/cr m < 1, these columns list NaN (see Section[44]l; 

- Columns 16-18: if, for the given value of the confidence level CQauss, the argument of the square root in Equation {20\ is negative, these 
columns list NaN (see Section [5~T| (; 

- Columns 16-18 and 20-22: if the radial- velocity precision (standard error) to reach a certain confidence level exceeds 9999.99 km s~' , a value 
of 9999.99 km s-' is listed. 
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